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Quantum mechanics is already 100 years old, but remains 
alive and full of challenging open problems. On one hand, 
the problems encountered at the frontiers of modern theo- 
retical physics like Quantum Gravity, String Theories, etc. 
concern Quantum Theory, and are at the same time related 
to open problems of modern mathematics. But even within 
non-relativistic quantum mechanics itself there are fundamen- 
tal unresolved problems that can be formulated in elementary 
terms. These problems are also related to challenging open 
questions of modern mathematics; linear algebra and func- 
tional analysis in particular. Two of these problems will be 
discussed in this article: a) the separability problem, i.e. the 
question when the state of a composite quantum system does 
not contain any quantum correlations or entanglement and b) 
the distillability problem, i.e. the question when the state of 
a composite quantum system can be transformed to an en- 
tangled pure state using local operations (local refers here to 
component subsystems of a given system). 

Although many results concerning the above mentioned 
problems have been obtained (in particular in the last few 
years in the framework of Quantum Information Theory), 
both problems remain until now essentially open. We will 
present a primer on the current state of knowledge concern- 
ing these problems, and discuss the relation of these problems 
to one of the most challenging questions of linear algebra: the 
classification and characterization of positive operator maps. 

I. INTRODUCTION 

Quantum Mechanics celebrates in this year its first cen- 
tury of life. In October 1900, Max Planck presented to 
the "Deutsche Physikalische Gesellschaft" his seminal pa- 
pers: "Uber eine Verbesserung der Wienschen Spektral- 
gleichung" and "Zur Theorie des Gesetzes der Energie- 
verteilung im Normalspektrum" Q]. 

While there arc no doubts about the success of Quan- 
tum Mechanics in explaining - beautifully - many of the 
problems concerning a variety of physical topics, it is 
worth stressing that Quantum Theory is by no means 
passee. It is nowadays still an open theory with several 
challenges which compel both Physics and Mathematics. 
Just a few years ago, it was common belief that the "really 
big" unsolved problems of theoretical physics pertained 
only to the domain of Quantum Gravity, that is the con- 
junction of Quantum Field Theory and General Rela- 



tivity. Quantum Gravity and String Theory are closely 
related to unsolved challenges of modern mathematics, 
in particular concerning algebraic topology and algebraic 
geometry. However, to the surprise of many, the emerg- 
ing field of Quantum Information Theory has shown 
that even in the "simple" non-relativistic Quantum The- 
ory there exist still fundamental open problems. The 
characterization of entanglement, and more specifically, 
the characterization of separability and distillability of 
quantum states are among these. Again, they are directly 
linked to unsolved challenges of mathematics concerning 
linear algebra and geometry, functional analysis and, in 
particular, the theory of C*-algebras 

In this paper we present some of these new open prob- 
lems and report on the recent progress concerning them. 
The paper does not intend to be a review article on the 
subject of quantum entanglement, but rather an intro- 
duction - a primer - on the subject|^ It contains nev- 
ertheless some new results: we present here two novel 
separability checks, and some new results concerning dis- 
tillability of density matrices that possess a non-positive 
partial transpose. The presentation of these new results 
and their proofs require the introduction of some tech- 
nical formalism. Therefore, they have been included as 
Appendices. 

The paper is organized as follows: First, in Section II 
we explain what the entanglement problem means. In 
Section III we discuss the problem of separability, that is 
how to define and discriminate those states that contain 
only classical correlations and no quantum correlations. 
In Section IV we focus on the problem of distillability, 
that is, the possibility- by using local operations and 
classical communication only- to "distill" from a given 
ensemble of copies of a given mixed state a maximally 
entangled pure state. In that section we also present the 
current state-of-the-art. Finally, our summary remarks 
are contained in Section V. 



^This paper has been presented by M. Lewenstein at the 
conference "Quantum Optics", Kiihtai, January 2000. 



1 



II. THE ENTANGLEMENT PROBLEM 

In order to explain what the entanglement problem 
means we first have to specify that in the following we 
will consider composite quantum systems Q. Physical 
states of such systems are in general mixed and can be 
represented by density matrices, i.e. hermitian, positive 
definite linear operators of trace one, acting in the Hilbert 
space H = TCa ® TCb ® • • ., which is a tensor product of 
Hilbert spaces corresponding to subsystems A,B,... of 
the considered system. 

Given a quantum state p G 7i, an apparently innocent 
question as does this state contain quantum correlations ? 
will in general be very hard (if not impossible!) to answer. 
First of all, what does it mean that a quantum state does 
or does not contain quantum correlations? The answer 
seems to be straightforward: a system contains quan- 
tum correlations if the observables of the different sub- 
systems are correlated, and their correlations cannot be 
reproduced by any means classically. That implies that 
some form of non-locality jj] is necessary in order to ac- 
count for such correlations. For pure states - described 
by projections on a single vector acting on the Hilbert 
space H of the composite system - it is relatively easy 
to check if the correlations that they contain are classi- 
cal, or not. For instance, it is enough to check if some 
kind of Bell inequality is violated to assert that the 
state contains quantum correlations. In fact, there are 
many different "entanglement"- criteria and all of them 
reveal equivalent forms of the non-local character of the 
entangled pure states. For example, the demonstration 
that no local hidden variable (LHV) can account for the 
correlations between the observables in each subsystem 
is an equivalent definition of non- locality |^ . 

We know nowadays that these equivalences may fade 
away when one deals with mixed states. Contrary to a 
pure state, a mixed state can be prepared in many dif- 
ferent ways. The fact that we cannot trace back how it 
was prepared prevents us from extracting all the infor- 
mation contained in the state. As a consequence, we lack 
(nowadays) general "entanglement" criteria that allow us 
to check if the correlations present in the system are gen- 
uinely quantum, or not. Despite the fact that many en- 
tanglement measures have been introduced, we do not 
know a "canonical" way of quantifying the entanglement 
. Furthermore, different manifestations of non-locality 
are known to be not equivalent. For instance, Werner 
introduced a family of mixed states that do not violate 
Bell-type inequalities (they admit a local hidden vari- 
able model), but nevertheless are non-local. The ques- 
tion whether there exists a violation of Bell inequalities 
in the, so-called, strong sense (where the observables take 
"unphysical" values) for all PPT entangled states (which 
we define below) remains open . 

Therefore, the entanglement problem can be outlined 



as: What does it mean that a given (mixed) state p con- 
tains or does not contain quantum correlations? 

III. THE SEPARABILITY PROBLEM 

An essential step forward to understand what does en- 
tanglement mean is to discriminate first the states that 
contain classical correlations only (or no correlations at 
all). These states are termed separable states, and their 
mathematical characterization has been formulated by 
Werner We shall restrict ourselves here to the most 
simple composite systems: bipartite systems (with two 
subsystems traditionally denoted as Alice and Bob) of 
finite, but otherwise arbitrary dimensions. The states 
of bipartite systems are described by positive definite 
hermitian density matrices (with normalized trace) p, 
i.e. p > 0, p = p^ and Trp = 1. The density ma- 
trices act on the Hilbert space of the composite system 
Ti. = 7i_4 (g) Tis- Without loosing generality we will as- 
sume that dim Ha — M > 2 and dim Hb — N > M . 

The most simple examples of separable states are just 
product states, i.e. p = <S) {p^ acts on l-LA^ and p^ 
acts on TYb). These states contain no correlations what- 
soever. A straightforward extension of product states are 
the states that contain only classical correlations. Werner 
Q provided us with the following operational definition 
of separability: 

Definition 1 A given state p is separable if and only if 

k 

1=1 

where X^iPi — 1; and pi > 0. 

The above expression means that p can be written as a 
convex combination of product states. Equation (1) has 
a clear physical meaning. The state p can be prepared 
by Alice and Bob by means of local operations (unitary 
operations, measurements, etc.) and classical communi- 
cation (LOCC). If p is separable the system does not con- 
tain quantum correlations. In spite of the definition, the 
characterization of such states is a rather arduous task. 
This is so among other facts because, in general, even 
for a given separable generic matrix we do not have an 
algorithm to decompose it according to Eq. (^. Thus, 
the separability problem, perhaps even more basic and 
fundamental than the entanglement problem can be for- 
mulated as: Given a composite quantum state described 
by p, is it separable or not? 

Before proceeding further we introduce here the defi- 
nitions that we shall use throughout the paper. Given a 
density matrix p, we denote by K{p), R{p), and r(p) the 
kernel, the range and the rank of the matrix p defined as: 

Definition 2 Kernel K{p) = {| 0) : p| = 0}. 
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Definition 3 Range R{p) = {| 0) : 3| V') : | </>) = p| -0)}. 

Definition 4 Rankr{p) — dim R{p) = N M — A\m K {p) . 

Let us introduce also the operation of "partial transpo- 
sition" that will be used throughout the paper and is 
defined as: 

Definition 5 The partial transpose of p means the trans- 
pose only with respect to one of the subsystems. If we ex- 
press p in Alice's and Bob's orthonormal product basis: 

M N 

i,j k,l 
M N 

= J2Y.(hk\p\j,l)\i}Aj\^\k)B{ll (2) 

i.j k,l 

then the partial transposition with respect to Alice is ex- 
pressed as: 

M N 
i,j k,l 

Note that p'^^ is basis-dependent, but its spectrum is not. 
The partial transpose p^-^ may be > 0, but does not have 
to be! As (p"^"*)"^^ = P"^, and as > always holds, 
positivity of p"^^ implies positivity of p"^^ and vice versa. 

A major step in the characterization of the separa- 
ble states was done by Peres Q and the Horodecki 
family [p^ . Peres provided a "userfriendly" and very 
powerful necessary condition for separability. Later on, 
Horodecki's demonstrated that this condition is also suf- 
ficient for composite Hilbert spaces of dimension 2x2 
and 2x3. Their results are enclosed in the following two 
theorems: 

Theorem 1 If p is separable then p^^ > 0. 

A matrix that verifies the above theorem is termed 
"PPT" for positive partial transpose. Notice that being 
a PPT state is a necessary condition for separability. 

Theorem 2 // p'^^ > in spaces of dimensions 2 x 2 or 
2x3 then p is separable. 

In general, there exist PPT states p (i.e. states with 
p^^ > 0) which are not separable in M x N spaces 
(Af = 2, > 4 or M > 3) 0. The PPT entangled 
states have been termed "bound entangled states" to dis- 
tinguish them from the "free entangled states" . This 
latter names are associated with the distillability prop- 
erty, which we will discuss in the later sections of this 
paper. "Bound entangled states" are entangled, how- 
ever, no matter how many copies of them we have, these 



states cannot be "distilled" via local operations and clas- 
sical communication to the form of a pure entangled state 
]l2t . We encounter thus new problems such as: How can 
one distinguish a separable state p from a PPT state p? 
Are all non-PPT states (NPPT states) "free entangled" 
i.e. distillable? But before trying to answer these ques- 
tions (that will bring us directly to the problem of dis- 
tillability) , we will first present a physical explanation of 
what separability means, and discuss shortly the recent 
progress concerning the quest for separability criteria. 

A. Physical interpretation of separability 

Let us now interpret the condition of positive partial 
transposition from a physical point of view. We start 
by considering symmetry transformations in the Hilbert 
space of each subsystem. Wigner's theorem tells us 
that every symmetry transformation is necessarily imple- 
mented by a unitary ([/) or anti-unitary (A) matrix. The 
tensor product of a unitary and an anti-unitary trans- 
formation Ua <E) Ab (or A A <8) Ub) results in a trans- 
formation which is neither unitary, nor anti-unitary in 
Ji = Ha ® Hb, and whose action on a general ket of 
the composite system [ip) € 7i, furthermore, cannot be 
properly defined. However, its action on a product ket 
I e, /) EE I e) ® I /), (where \ e) ^Ha and | /) G Hb) is, 
apart from a phase ambiguity, well defined. Thus, the ac- 
tion of a combined transformation of the type Ua (X" Ab 
on projectors corresponding to pure product states is well 
defined without any ambiguity. As a separable state p^ 
can always be rewritten as a statistical mixture of prod- 
uct vectors (see Def. 1) it is clear that under the com- 
bined transformation Ua <^ Ab (or A a <8) Ub), Ps trans- 
forms into: 

k 

p^^/'l = E?''(i^i)(^ii®i/^)</^i) w 

i=l 

where | e^) = UaIc,) € Ua ; \ f[) = e Ub- 

Therefore, p^ describes also a physical state so that p^ 
is a positive definite hermitian matrix (with normal- 
ized trace). This is what characterizes separable states: 
that any local symmetry transformation, which obviously 
transforms local (in this context local refers to each of the 
subsystems) physical states into local physical states, also 
transforms the composite global state into another phys- 
ical state. 

There exists only one independent anti-unitary sym- 
metry [ p^ , and its physical meaning is well known: time 
reversal. Any other anti-unitary transformation can be 
expressed in terms of time reversal (as the product of 
a unitary matrix times time reversal). Thus separabil- 
ity of composite systems implies the lack of correlation 
between the time arrows of their subsystem. In other 
words: given a separable composite state, reversing time 
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in one of its subsystems leads again to a physical state 
00 • 



B. Quest for separability criteria and checks 

In the recent years there has been a growing effort in 
searching for necessary and sufficient separability crite- 
ria and checks. Several necessary conditions for separa- 
bility are known: Werner has derived a condition based 
on the analysis of local hidden variables (LHV) models 
and the mean value of the, so-called, flipping operator 
the Horodecki's have proposed a necessary criterion 
based on the so-called a-entropy inequalities jl^, etc... 
Quite recently, a general and sufficient condition for sepa- 
rability was discovered by the Horodecki family in terms 
of positive maps. A map is defined positive if it maps 
positive operators into positive operators. The condition 
found by the Horodecki's, states that p is separable iff 
the tensor product of any positive map acting on one 
subsystem A and the identity acting on other subsystem 
B maps p into a non-negative operator. This definition, 
however, involves the characterization of the set of all 
positive maps which is per se a major task. Later on the 
reduction criterion of separability was introduced 0,|T||: 

Criterium 1 If p is separable then the map (Ta <Xi 1) : 
p (1a ® TrAp) — p must be positive. 

Violation of this criterion is sufficient for entanglement 
to be free. Following the reduction criterion, a simple 
and still quite powerful sufficient condition for distilla- 
bility was provided in Rcf. ||T^, where P. Horodecki et 
al. showed that if the rank of at least one of the reduced 
density matrices pA = Tr^p and pb — Tr^p exceeds the 
rank of p, then p is distillable, ergo is non-separable and 
NPPT. In particular, it was concluded in Ref. that 
there is no bound entanglement of rank 2. 

Sufhcient conditions for separability are also known. 
In Ref. it was proven that any state close enough to 
the completely random state tt = I/NM is separable. In 
|pO| there were also given the first quantitative bounds for 
the radius of the ball surrounding tt that does not contain 
any entangled state. Much better bounds were found in 
the following works [pl| , where it was proven that a full 
rank mixed state is separable provided that its smallest 
eigenvalue is greater or equal to (2 -I- MN)~^. 

In Ref. , in which the first explicit examples of en- 
tangled states with PPT property were provided, another 
necessary criterion of separability was formulated. Ac- 
cording to this criterion: 

Criterium 2 // the state p acting on a finite dimen- 
sional Hilbert space is separable then there must exist a 
set of product vectors {jci, /i)} that spans the range R{p) 



such that the set of partially complex conjugated product 
states {|e*,/i)} spans the range of p"^^ . 

The analysis of the range of the density matrices, initi- 
ated by P. Horodecki, turned out to be very fruitful, lead- 
ing, in particular, to the algorithm of optimal decompo- 
sition of mixed states into the separable and inseparable 
part 1^,^ , and to systematic methods of constructing 
examples of PPT entangled states with no product vec- 
tors in their range, using either so-called unextendible 
product bases (UPB's) or the method described 

in (26). 

In the Appendices A and B we present two novel sep- 
arability criteria and checks. One provides a necessary 
condition for separability, or rather a sufficient condition 
for entanglement. It detects the non-separability of the 
UPB states. The other criterion, or rather separability 
check, detects all separable states that are the convex 
combinations of two product states. 



C. Recent progress in the separability problem 

Despite many efforts and seminal results obtained in 
the recent years, the problem of separability remains es- 
sentially open. Recently, a considerable progress of in 
the study of PPT entangled states has been made |^,^ . 
The results obtained allow us to hope to develop a sys- 
tematic way of constructing optimal criteria for separa- 
bility in arbitrary Hilbert spaces |^ . 

Our method employs the idea of "subtracting projec- 
tors on product vectors" [ p2p^ ] : if there exists a product 
vector |e, /) e R{p) such that |e*,/) € R{p'^^), the pro- 
jector onto this vector (multiplied by some A > 0) can 
be subtracted from p, such that the remainder is positive 
definite and PPT. Our results can be divided into three 
groups. 

First, we have studied and found separability crite- 
ria for density matrices of sufhciently low dimensional 
rank. Also constructive algorithms to decompose opti- 
mally (with the smallest possible number of terms) the 
separable matrix according to Eq. (1) - i.e. in product 
states - have been provided for low rank matrices. For 
the general case of composite systems with the Hilbert 
space (g) Ti^ (M < N) our findings are essentially 
contained in the following two theorems: 

Theorem 3 If p is PPT such that r(p) = k < N and p 
cannot be embedded in a 2 x (k — 1) dimensional space 
then p is separable. 

In particular when p has rank N there exist typically 
exactly N product vectors jci, fi) in the range of p such 
that |e*, fi) £ R{p^^); p is then a convex combination of 
projections onto these vectors. 
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Theorem 4 // r{p) + rip^^ ) < 2NM - N-M + 2 then 
typically there exists a finite number of product vectors 
I e, /) e R{p) such that \e*J)e R{p^^). 

These product vectors are the only possible candidates 
to appear in the decomposition of Eq. (1). Finding them 
requires solving a system of polynomial equations. After 
these equations are solved, one can check whether p has 
the decomposition (1). The problem is infinitely easier 
than the original one since we know now all possible pro- 
jectors that can be used, and we know that their number 
is finite. In fact, checking in such a situation whether 
p is separable or not can be done in a finite number of 
computational steps! 

Second, we have studied the structure and generic form 
of low rank PPT entangled matrices. To study low rank 
PPT entangled matrices has a twofold purpose. On one 
hand, the complexity of the problem is reduced, and 
therefore it is possible to find separability criteria. But, 
perhaps the most important is the fact that, given a den- 
sity matrix p, one can always decompose it as : 

p = Kpsep + (1 - A)(5p , (5) 
p^-^Kp^\,.^ + {l-K)5p^- , (6) 

where ps^p — X^i ^iPi is a separable state. Pi are pro- 
jectors onto product states and A is maximal. All the 
information concerning entanglement is then contained 
in the remainder 5p {5p'^^) which has low rank and can 
be termed as a "pure" PPT entangled state, or "edge" 
PPT entangled state. This state has a property that no 
projection onto the product state can be subtracted from 
it, keeping the rest positive definite and PPT. Formally, 
there exist no product vectors |e,/) G R{Sp) such that 
|e*,/) e R{6p^^). The "edge" states violate in the ex- 
tremal sense the Criterion of Rcf. . The problem of 
the separability reduces now to the problem of separa- 
bility of the "edge" states, and to the question whether 
a given mixture of an "edge" and a separable state is 
separable or not. 

In other words, any matrix p can be decomposed in 
a separable part (that contains product vectors in the 
range) and a remainder, which is the "edge" state. Sp 
and Sp^^ are low rank matrices that contain all the in- 
formation related to entanglement. Obviously, knowing 
the structure of those matrices is therefore of capital im- 
portance. 

Finally, let us mention a different approach to the en- 
tanglement problem, that is based on the so-called en- 
tanglement witnesses. An entanglement witness is an 
observable E that reveals the entanglement of an entan- 
gled density matrix p. B. Terhal introduced entan- 
glement witnesses through the following theorem: 

Tiieorem 5 If p is entangled then there exists an entan- 
glement witness E such that 



T^r{Ep,ep) > 0, (7) 
Tt{Ep) < (8) 

for all separable matrices Psep- 

Entanglement witnesses represent - in some sense - a 
kind of Bell inequality which is violated by the entan- 
gled state p. Each entanglement witness i? on an M x 
space defines a positive map £ that transforms positive 
operators on an M or Af-dimensional Hilbcrt space into 
positive operators on an M or A^-dimcnsional space . 
The maps corresponding to entanglement witnesses are 
positive, but not completely positive, and in particular 
their extension to M x A^ spaces allows to "detect" the 
entanglement of p. The maps corresponding to entangle- 
ment witnesses for PPT states are, moreover, non decom- 
posable: they cannot be represented as a combination of 
completely positive maps and partial transposition. 

For every "edge" state it is possible to construct an 
entanglement witness ||2^ as: 

E = PK(,) + {P^^pT./^-^, (9) 

where Px{p), Pj^(^pTA^ a-^e projections onto the kernel of 

p and the kernel of p"^-*, and < e = min(e, f\PK{p) + 
(P^^^Ta^"^"* |e, /), where the minimum is taken over all 
possible product vectors. 

We have not only been able to find entanglement wit- 
nesses, and the corresponding non-decomposable posi- 
tive maps for arbitrary "pure" or "edge" PPT states, 
but also to optimize them in a certain sense |^^. Opti- 
mized entanglement witnesses detect significantly more 
entangled PPT states than the non-optimized ones. 

We hope very much that these studies will allow us to 
characterize extremal points in the convex set of PPT en- 
tangled matrices, and then to characterize the extremal 
points in the convex set of positive maps ||3C| ] . If this pro- 
gram is realized, the separability problem will be solved. 
So far, however, only the first steps have been done and 
the problem remains open and challenging. 

IV. THE DISTILLABILITY PROBLEM 

On having said that, we shall attack now the related 
problem of the distillability of mixed quantum states. 
For many applications in quantum information process- 
ing and communication one needs a maximally en- 
tangled state, that is, a state which in Mx A^ dimensional 
space can be brought by a local change of basis to the 
form 

1 

= I z,z), (10) 

^ i—l 

which is shared between two parties. 
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Although in principle one can create pure and max- 
imally entangled states, in realistic situations any pure 
state will evolve to a mixed state due to its interac- 
tion with the environment. A standard example con- 
cerns a situation when two entangled particles (photons, 
atoms,...) representing the two subsystems are sent from 
the source to the two involved parties, Alice and Bob, 
through noisy channels. In order to overcome the noise 
created during the transmission, the idea of distillation 
and purification, i.e. enhancement of the given non- 
maximal mixed entanglement by local operations and 
classical communication (LOCC) was proposed by Ben- 
nett et al. [Q, Deutsch el al Q and Gisin Again, 
for Hilbert spaces of composite systems of dimension 
lower or equal to 6, any mixed entangled state can always 
be distilled to its pure form. Since for such systems en- 
tanglement is equivalent to the NPPT property, we con- 
clude that for 2x2, and for 2x3 systems, all NPPT states 
are distillable |Q. However, it was shown by Horodecki 
family that in higher dimensions there exist states 
(namely PPT entangled states), termed as bound entan- 
gled states, which cannot be distilled, in contraposition 
to free entangled states |3^. The distillability problem 
can be formulated as: Given a density matrix p, is it or 
is it not distillable? 

Let us now define the distillability property, first on an 
intuitive, and then on a more formal basis. 

Definition 6 p is distillabe if by performing LOCC on 
some number K of copies p, Alice and Bob can distill a 
state arbitrary close to \ '^max), *-e- 

p®....(^p >\'i!rnax) {"UmaxV (H) 

The above definition is not very precise - it requires to 
specify what the LOCC can do with K copies of p, and 
does not give any practical advice about how to answer 
the question of distillability. Fortunately we can use the 
theorem of Ref. , which states that instead of studying 
the whole set of possible LOCC, in order to determine 
the distillability of a given density matrix it is sufficient 
to study projections on a 2 x 2 dimensional subspace of 
the Hilbert space in which p acts. The theorem is very 
useful since it reduces the problem of distillability to a 
very precisely stated mathematical question; in fact from 
now on we will use it as a definition of distillability. 

Theorem 6 p is distillable iff there exists a number of 
copies K , and a projector P2x2 onto a 2 x 2-dimensional 
space spanned by: 

I e,)e HA^-^nA , * = 1, 2 , (12) 

^ v 

K —times 

\fr)€nB(E>...(E>nB , 1 = 1,2, (13) 

K — times 



such that the projection 

<T = P2X2P^'^P2X2, (14) 

is NPPT (i.e. is distillable). 

An alternative way of formulating the above theorem is 
the following: p is distillable iff there exists a state | ip) 
from a 2 x 2-dimensional subspace, 

|V') = a|el)|/i)+6|e5)|/2) , (15) 

such that {ip |(p'^'*)^^| -0) < for some K. 

The idea of the proof of the above theorem is the fol- 
lowing: if p is distillable it means that one can produce a 
maximally entangled state, and it is then easy to project 
(using local projections) that state onto a pure state in a 
2 X 2-dimensional subspace. On the other hand, if there 
is a as in equation (|l5|), such that \p'^'^\'ip) < 0, 
then one can first project onto the 2x2 subspace to 
which I ip) belongs. This is a 2 x 2 subspace in which 
the projected matrix is NPPT, ergo it is distillable. We 
can then distill several maximally entangled state in this 
subspace, rotate them unitarily and locally, and combine 
to a maximally entangled state in the whole space. 

Let us now ask what does the criterion of partial 
transposition,- which plays an important role in the sep- 
arability problem as we have seen before-, tell us about 
the distillability problem? 

Theorem 7 If p"^^ > (p is PPT) then p is not distill- 
able 

Theorem 8 //p^-* (pis NPPT) in dimensions 2x2, 
2x3 then p is distillable ^d^. 

The later holds also for 2 x systems, see js^ . 

A. Recent progress in the distillability problem 

At the end of the last section we have seen that every 
density matrix with a positive partial transpose cannot 
be distilled, and that for low dimensions the converse 
is true. In this section we want to discuss the conjec- 
ture that in higher dimensions there are states with a 
non-positive partial transpose, which are however non- 
distillable [^,^. In other words, non-positivity of the 
partial transpose seems to be a necessary, but not a suf- 
ficient condition for distillability. 

The states which are believed to be non-distillable , 
belong to a one-parameter family, which lives in dimen- 
sion N X N: 

p{a) = —^{Ps + aPA) with a > , (16) 
mia] 
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where Ps and Pa denote projectors onto the symmetric 
and antisymmetric subspace, respectively, and m is some 
normalization. This family is generic in the sense that 
every density matrix in N x N can be depolarised locally 
to a state from our family. In this sense these states are 
nothing more but Werner states ||^ defined for N x N 
systems. 

The partial transpose of p{a) is given by 



1 



(1 - pP) 



(17) 



where n(/3) is the normalization, P is the projector onto a 
maximally entangled state l^'maa; 
the relation between a and /3 is: 



Eti \^^)/VN, and 



f3- 



N{a-1) 
a + 1 



with - N < ^ < N 



(18) 



Note that for /3 < 1 the matrix p^^ is positive definite , 
i.e. p is PPT and thus is not distillable. For (3 > 1 one 
finds that p^-* is not positive definite, i.e. p is NPPT and 
the question of distillability is open. 

The nice thing about the considered family of states is 
that if we show that: 

(i) p is distillable for all /3 > 1, then all p with NPPT are 
distillable, because all p can be reduced to the "canoni- 
cal" form (|l6|). 

(ii) if there exists P < N, such that p is NPTT and is 
not distillable, then not all p with NPPT are distillable. 
In other words there exist undistillable p's with NPPT. 
Both alternatives (if proven) would be an extremely im- 
portant result. At the moment it seems that the second 
alternative is true |^,^ , but strictly speaking the prob- 
lem is open. We shall see below that p is distillable for 
P > 3/2, so that in fact interesting region of the param- 
eter P lies between 1 and 3/2. 

Before presenting some partial results on the way to 
the complete proof, which is yet unknown, let us intro- 
duce the concept of if -distillability, by which we define 
distillability with respect to K copies. 

Definition 7 p is K -distillable iff there 3\ ip) as given in 
Eq. ( |7^ j such that 



(^|(p^-)«^|V) <o. 



(19) 



p is not 2- distillable for 1 < /? < |; 



• p is K-undistillable for 1 < P < Pk , where the best 
bound for Pk obtained so far is: 



K 



1 



1 



3^/3/^1/3 ■ 



(20) 



As we see, for every K there exist a region of p (see 
Fig.l) in which p is not distillable. Unfortunately, this 
region shrinks, however to a point, as K goes to infinity. 
If we have an arbitrary number of copies of p we cannot 
say whether we will be able or not to distill it. The 
proof of the above theorem is technical, but essentially 
simple. Below we sketch the proof concerning 1- and 
2-distillability 



B. 1-Distillability 

Let Q = 1 — P be the projector complementary to 
P, where P is the projector onto maximally entangled 
states. For one copy of the matrix p we have that for any 
I ?/;) given by Eq. (|TF 



{^\Q-pP\^) = {i;\l-{l+P)P\^)> 
i-{l + P)\>Q 



(21) 



for /3 < |. The last inequality follows from the fact that 
the projection |('0|^'maa;)P of an arbitrary vector living 
in a 2 X 2 subspace onto a maximally entangled vector 
in 3 X 3 space must not be greater than 2/3. The above 
considerations imply also that for P > 3/2 the matrix 
is 1-distillable, because there exists a vector in a 2 x 2 
subspace for which KV'l^maa;)^ — 2/3. This proves that 
p is distillable for P > 3/2. 



C. 2~Distillability 

To show that 2 copies of p cannot be distilled for some 
interval of P, we first observe that Q is separable, i.e. 



Q ^^Pi\ei, fi){et, fi 



(22) 



The basic theorem that we have proven so far determined 
the region of the parameter P for which the matrix p 
of Eq. ( pj| ) is not K-distillable. Our results hold for 
arbitrary N x N systems, but here we specify them to 
case of two qutrit systems (3x3). 



Theorem 9 Let p ^ Ps 

m 3 X 3 space, then: 



uPa, i-e. p'^^ ~ 1 — pp act 



• p is not 1-distillable for 1 < P < 



Denoting by Tr^ the trace over the i-ih copy, and by Qi 
- the projector Q for the i-th copy we get 

R 

Tr((V' IQil V^)) = Y^pM /'I ^) (23) 

i=l 

where = ("0164,/^) is a vector in the second copy 
space with 2 Schmidt coefficients, i.e. of the form (p^). 
Similar result holds for the second copy. Then, using the 
results for 1-distillability we obtain 



7 



{^\Q®iQ~^P)\i^)>0, 
{^\iQ~lp)(^Q\i')>0- 



(24) 
(25) 



- separable - 



Now, by adding the above results and dividing by two we 
obtain 

(V- |(g - ip) ® (Q - ip)| V) > I^P ® P)\ i') > 0, 

(26) 

so that we see that p is not 2-distillable for P < 5/4. 



d+l 

d-1 



- pos. part, transp. 



- noivpos. part, transp. 



FIG. 1. Separability and distillability properties of p in 
arbitrary d x d dimensional space. In the region between 
separable and distillable p is NPPT and most presumably 
non-distillable. 



D. Distillability in general 



APPENDIX A: SUFFICIENT CRITERION FOR 
INSEPARABILITY 



I we have performed extensive nu- 
and looked for the minimum of 



In the Ref. [ 
merical studies 
('0 |(p^'*(/3))®'^| V') over all possible | ip) of the form 1^ 
The numerical results indicate clearly that: 



• p is not 2-distillable for 1 < /3 < |; 
9 p is not 3-distillablc for 1 < /? < |. 



It is a challenging and open problem to understand these 
results. So far we have only achieved some progress in 
the problem of 2-distillability. We have proven that 
the states \4>) for which (V' |(p'^^(/3))®-^| V') = 0, and 
which (as we know) provide the global minimum of 
(-0 Kp^-^ (/3))®^| V') for /3 < 5/4, provide also a local min- 
imum of {^p\{p'^^{P))'^^\-)p) equal to zero for 3/2 > /3 > 
5/4. The proof of this fact is presented in the appendix 
C. There exists a well-based suspicion that p is not dis- 
tillable in the entire region ofl</3<3/2. 



V. CONCLUSIONS 

There is only one conclusion of this paper: Quantum 
Theory is an open and challenging area of physics. It of- 
fers still fundamental and fascinating problems that can 
be formulated at elementary level, and yet they are re- 
lated to challenges of the modern mathematics. Partic- 
ular examples of those are separability and distillability 
of composite quantum systems. 

This paper has been supported by SFB 407 
and Schwerpunkt "Quanteninformationsverarbeitung" of 
Deutsche Forschungsgemeinschaft, by the ESP PESC 
Programme on Quantum Information, and by the 1ST 
Programme EQUIP. 



Our sufficient criterion of inseparability (i.e. a neces- 
sary criterion for separability) is based on the following 
Lemma. Let us assume the product vectors |e,/) to be 
normalized. 



Lemma 1 // maa;{Tr(p2),Tr((p'^«)2),Tr(pp'^«)} 
maa;|e./) (e, /l/oje, /), then p is inseparable. 



> 



Let us prove that this statement is true in the case where 
{Tr(p2) = max{Tr(p2),Tr((p^«)2),Tr(p/9^«)}. All the 
other cases can be proved in the same way which implies 
that we can take the maximum of those three values. 
Proof: We define r = max|ej) (e, /|p|e, /) and the wit- 
ness, X = 1 — ip. Assuming that Tr(p) = 1 we have that 
Ti{Xp) < 0. It remains to prove that Tt{Xp) > for all 
p separable. We write p = J2i fi){^i^ fi\^ ^^'^ ob- 
serve that: rTr(Xp) = Tr[(max|ej) (e, /|p|e, /) - p)p] = 
J2i Ai(max|ej)(e,/|p|e,/) - {eiji\p\e,ji)) > 0. 

Unfortunately this criterion does not work for the 
Horodecki PPT states in 2 x 4 space ||l^]. It does, how- 
ever, detect the entanglement of the PPT states con- 
structed from UPB's [Q. In such case p — P/K, where 
P is a projector onto a space that does not contain any 
product vectors, K = r{P), and P — P^'^ . Our crite- 
rion gives Tr(p2) ^ Tr((p^s)2) = Tr(pp^«) = 1/if, while 
(e,/|p|e,/) = {eJ\P\eJ)/K <1/K. 



APPENDIX B: SEPARABILITY CHECK FOR 
BINARY MIXTURES OF PRODUCT STATES 

If a separable matrix is a mixture of two product states 
(here we call such matrix a binary mixture), then it is 
relatively easy to check separability. Assume p of the 
following form: 



K 



(Bl) 



i=l 
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where all pi > 0, and J^iPi = 1- Let and i/ be density 
matrices acting in Alice and Bob's space respectively. Let 
us define the matrix function 



K 



Mill, v) = Y, Mpf -^^)® {pf - 



(B2) 



Interestingly, M can be calculated without the explicite 



use of the representation (Bl), 

M(p, v)^p~^,(i^pB-pA®v + ^^v, (B3) 

where the reduced density matrices are pa,b — Trs,A/0- 
For K — 2 'we have an obvious Lemma: 

Lemma 2 If p is a mixture of two product states ( with 
Pi — p,P2 ^ I — p) then the equation M{p,v) — has at 
least two solutions p = pf, v ~ p2 , o-nd p — p2, v ^ pf . 

The "opposite" implication is also true. 

Lemma 3 // the equation M{p, v) = Q has solutions 
such that p > Q, p — p\ Tr/i — l,v>Q,v — v\ 
Tri/ = 1, then p is either a separable binary mixture, or 
a nonseparable binary pseudomixture. 



The proof follows directly from Eq. (B3). p can be writ- 
ten as 



p = p®{pB- {I- p)v) + [pA - pp) ® V, 



(B4) 



with some < p < 1. We immediately see that p is 
separable if there exists p € [0, 1] such that pB — {^~p^v > 
and PA — PP~> 0. This impHes that the ranges R{p) 
and R[v) must be included in the ranges of pA-, and ps, 
respectively; at the same time p must fulfill the conditions 
p\\p{pA)-^\\ < 1, and {I - p)\\pipB)-'\\ < 1, where ||.|| 
denotes the operator norm (for details see |p7||). 

Checking if the equation M{p, i/) = has solutions 
is very easy. We can use a product basis {Of ® 
0^}i=i_..._M2y"=i,...,Ar2 in the operator space. Such ba- 
sis can be chosen to be orthonormal and hermitian 
with respect to the trace scalar product. The equation 
M{p, i^) = projected onto the «j-th element of the basis 
reads: 



Pij - PiPBj - PAiVj + PtVj = 0. 



(B5) 



where pij = TipO^ 



Of, PM 



PBj 



TrpBOf, etc. ..We have thus N^AI^ such equations for 
+ N'^ real coefficients pi, Vj. The equations have a 
very simple structure and therefore it is easy to check: 
a) if they have a solution; b) if the resulting p and v are 
positive definite; c) if there exists p such that both terms 
on the RHS of Eq. (B4) are positive definite. 

The above formulated separability check can be easily 
generalized to systems of R parties and separable mix- 
tures of R product states for i? > 2. 



APPENDIX C: FINDING THE LOCAL MINIMUM 
FOR PROJECTING ONTO THE 
2-DIMENSIONAL SUBSPACE 

In this appendix we will find the states that lead to 
a local minimum for projection onto a two-dimensional 



subspace as in equation (19) in the case of 2 copies (in 
dimension 3x3), for the critical value of the parameter 
P = 3/2. 

We will proceed as follows: our problem will be for- 
mulated in terms of a function f{X,ip) that has to be 
minimized with respect to We will find a family of 
states ip for which / is shown to reach a local minimum, 
for a range of parameters A. 

We introduce a parameter A and will study the follow- 
ing function: 

/(A,V') = ^(^I(1 + |1-2A|)]1®1+ Jp®P 

-3(A-1®P+(1-A)P®1)|V) 

= (V'|P2(A)|V;) , (CI) 

where the last line defines /92(A), and P is the projector 
onto a maximally entangled state. Here A is a given fixed 
parameter with < A < 1. The case A = ^ corresponds 
to P2 being the partial transpose of two copies of 1 — 
3P/2. In the notation of Section IV this corresponds to 
the value of /3 = 3/2. We are looking for the minimum 
of / with respect to | . This state lives in the two- 
dimensional subspace and can be written in the Schmidt 
decomposition (cf. B) 



1^) =a|ei)A|/i)B 62)^1/2) 



B 



(C2) 



where the states | Ci), | fi) are normalised and (e| |e2) = 
= (/* 1/2) and I a |2 I 6 p = 1. For clarity we kept the 
indices A and B. 



Let us rewrite the terms in equation (CI), sorting them 
not pairwise, but with respect to Alice and Bob: 



1(g) 1 = 



P(X) P = - 



3 

E 

z,j,r,s — 1 
1 3 



P® 1 = - 



P = - 



3 

E 

i,j,r,s—l 
3 

E 



\ir)A\ir)B{j3\A{js\B, 



ir)A\is)B{jr U(js |b, 
ir)A\jr)B{is\A{js\B 



(C3) 



This notation fixes the basis in which we will also write 
I ij}). Indices i,j are used for the first pair and r, s for the 
second pair. 

The minimum of /(A 0, tjj) is found by requiring the 
two conditions 
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a) {i}\P®P\i}) = 0, 

b) (-0 |P ® Ij -0) maximal 



According to equations (C3) we reach a) only if the en- 
tries of either the first or the second pair are orthogonal 
to each other. We can reach b) if the entries in the first 
bits of Alice and Bob are identical in both terms of the 
Schmidt decomposition of | V") , and if their second bits 
are in a product state. The coefficient in equation ( |C2| ) 
is easily found to be | a | = l/\/2 for maximisation of b). 

We can fulfill both conditions a) and b) with the family 
of states that minimizes / for A = 0, denoted by | -0*): 



r) = ^(1 ^r)A\^s)B + (C4) 

we have 



with — = (r|s). Therefore for A 

found the global minimum of / to be 

min/(A = 0) = . 



(C5) 



Using the explicit structure of the states | ■0*) we find 
that 

(0*11 110*) = 1 , 
(0*|1(0P|^*) = , 

(V'*|P(^1|0*) =2/3 . (C6) 

When varying the parameter A we therefore find 

(V'*|P2(0< A< i)|V/) =0 . (C7) 

Similarly, for ^ < A < 1 the same line of argument holds 
when interchanging the role of first and second bits, lead- 
ing to a different minimizing family | 0*). At A = ^, the 
point which is symmetric with respect to interchanging 
the two pairs, both families lead to the expectation value 
zero. 

In the following we will show that | 0*) corresponds to 
a local minimum for < A < i. 



First, we show that the states given in equation (C4) 
form a compact set, by describing how to move through 
the whole family in infinitesimal steps: Looking at the 
first pair, we can either make the following change: 

\i) ^ xM) + Xk\k) with (k |i) = = (k |j) , (C8) 

or we can move to 

\j)^x,\j)+xi\l) with(l|j>=0=(l|i) , (C9) 

or we can change both \ i) and keeping — 0. 

Regarding the second pair, we can change 

I r) ^ Xr\ r) + Xj,\p) with (p |r) = = (p |s) , (CIO) 

or we can move to 

\s) ^ Xs\s) +xt\t) with (t|r) = = (t|s) . (Cll) 



In this way we can move within the family in infinitesimal 
steps, and there are no isolated points. 

Let us now move outside of our family by an infinites- 
imal amount. We will write down the most general path 
leading away from the family and then show that first or- 
der terms of the expectation value vanish, i.e. we have an 
extremum, and that the functional determinant of second 
order terms is positive, i.e. we have a local minimum. 

The most general infinitesimal step away from our fam- 
ily is given by 



I 0-* + 5) = ■^Wl + 5o\ ir)A\ {i + 5ik){s + 52r))B + 



VT^e^'^l (j + 5Mr + 5At))A\ (j + 5zm){s + Jer))^) 
with (k|i) =0= (l|j), 
(/|z)=0 = (t|r), 
{m\]) ^ Q = {i + 5ik\j + S^m) , (C12) 

so that the Schmidt terms are still orthogonal, and for 
each of the seven (5i ^ we leave the family. Note that 
we can always keep one state (in this case | ir)^) constant 
by using bilateral rotations. 

We now expand the expectation value 



(^2) = ^(0* + (5 1^2(0 < A < i)| 0* + 5) (CIS) 



in powers of (5^, and find that all terms linear in 5i are 
indeed vanishing. 

The second order terms can be written down explicitly. 
The diagonal ones are: 



0{5l) : 




A) > 






0{5l) : 




1 1 - 2A|) > 






0{6l) : 




|l-2A| + i- 


1) > 




0{5l) : 




1 1 - 2A|) > 






0{5l) : 


5-54(1 + 


|l-2A| + i(5t, 


- (A5fs - 


h 1 - A)) > 


0{5l) : 




1 1 - 2A|) > 






0{5l) : 




|l-2A| + i- 


1) > . 


(C14) 



Here 5ts is the Kronecker symbol. 

Nearly all off-diagonal terms of second order vanish, 
the only non-zero one is [note that 0{5i5j) — 0{5j5i)]: 



(C15) 



the 



This term is negative for the range < A < ^; 
corresponding 2x2 determinant, however, is det(2, 6) = 
SA^ - 1 A-l- 1 > for < A < 1/2. Thus we have found the 
second derivative to be positive, and therefore our family 
I 0*) corresponds to a local minimum for < A < i. In 



particular that is the case for A 



i.e. 13 = 3/2 in the 
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notation of Section [V. It is easy to see that this must 
also be the case for all 1 > /3 < 3/2. 

For a complete proof that two copies are not distillable, 
however, it remains to be shown that this minimum is a 
global minimum for /3 = 3/2. 
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